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Abstract. Probing the fluid dynamics of thin films is an excellent tool to study the 
solid/liquid boundary condition. There is no need for external stimulation or pumping 
of the liquid due to the fact that the dewetting process, an internal mechanism, acts 
as a driving force for liquid flow. Viscous dissipation within the liquid and slippage 
balance interfacial forces. Thereby, friction at the solid/liquid interface plays a key role 
towards the flow dynamics of the liquid. Probing the temporal and spatial evolution 
of growing holes or retracting straight fronts gives, in combination with theoretical 
models, information of the liquid flow field and especially the boundary condition at the 
interface. We review the basic models and experimental results obtained during the last 
years with exclusive regard to polymers as ideal model liquids for fluid flow. Moreover, 
concepts that aim on explaining slippage on the molecular scale are summarized and 
discussed. 
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1. Introduction 

Understanding liquid flow in conflned geometries plays a huge role in the fleld of micro- 
and nanofluidics . Nowadays, microfluidic or so-called lab-on-chip devices are utilized 
in a wide range of applications. Pure chemical reactions as well as biological analysis 
performed on such a microfluidic chip allow a high performance while solely small 
amounts of chemicals are needed. Thereby, analogies to electronic large-scale integrated 
circuits are evident. Thorsen et al fabricated a microfluidic chip with a high density of 
micromechanical valves and hundreds of individually addressable chambers [2]. Recent 
developments tend to avoid huge external features such as pumps to control the flow by 
designing analogues to capacitors, resistors or diodes that are capable to control currents 
in electronic circuits [3j. 

By reducing the spacial dimensions of liquid volume in conflned geometries, slippage 
can have a huge impact on flow dynamics. Especially the problem of driving small 
amounts of liquid volume through narrow channels has drawn the attention of many 
researchers on slip effects at the solid/liquid interface. The aim is to reduce the pressure 
that is needed to induce and to maintain the flow. Hence, the hquid throughput is 
increased and, what is important in case of polydisperse liquids or mixtures, a low 
dispersity according to lower velocity gradients perpendicular to the flow direction is 
generated. Conflned geometries are realized in various types of experiments: The physics 
and chemistry of the imbibition of liquids by porous media is of fundamental interest 
and enormous technological relevance Pj. Channel-like three-dimensional structures 
can be used to artiflcially model the situation of fluids in conflnement. Moreover, the 
small gap between a colloidal probe and a surface flUed with a liquid (e.g realized in 
surface force apparatus or colloidal probe atomic force microscopy experiments) is a 
common tool to study liquid flow properties. Besides the aforementioned experimental 
systems, knowledge in preparation of thin polymer fllms has been extensively gained due 
to its enormous relevance in coating and semiconductor processing technology. Such a 
homogeneous nanometric polymer fllm supported by a very smooth substrate, as for 
example a piece of a silicon (Si) wafer, exhibits two relevant interfaces, the liquid/air 
and the substrate/liquid interface. Si wafers are often used according to their very low 
roughness and controllable oxide layer thickness. Yet, also highly viscous and elastically 
deformable supports, such as e.g. polydimethylsiloxane (PDMS) layers, are versatile 
substrates. 

The stability of a thin fllm is governed by the effective interface potential as 
function of fllm thickness h. In case of dielectric systems, (/>(/i) is composed of an 
attractive van-der-Waals part and a repulsive part [3, El [3, IH] . For the description of 
van-der-Waals forces of a composite substrate, the layer thicknesses and their respective 
polarisation properties have to be taken into account [9]. Thereby, three major situations 
of a thin liquid fllm have to be distinguished: stable^ unstable and metastahle fllms. As 
illustrated in Fig. [U by the typical curves of 0(/i), a stable liquid fllm is obtained if 
the effective interface potential is positive and monotonically decaying (cf. curve (1) in 
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Figure 1. Different shapes of the effective interface potential (j){h) associated with 
different wetting conditions. Curve (1) characterizes a stable liquid film. Curve (2) 
represents a metastable, curve (3) and (4) an unstable situation. 



Fig. [T]). The equilibrium film thickness h^g is infinite and the liquid perfectly wets the 
substrate. In case of a global minimum of 4>{h), curves (2) and (3) in Fig. [H the system 
can minimize its energy and a finite value for hgg results. The metastable situation is 
furthermore characterized by a potential barrier that the system has to overcome to 
reduce its potential energy (cf. curve (2) in Fig. [1]). Curve (3) and (4) characterize 
unstable conditions, since every slight fiuctuation in film height will drive the system 
towards the global minimum. The wettability of the substrate by the liquid is correlated 
to the depth of the minimum at (pijieq). The deeper the global minimum of 0, the larger 
is the equilibrium contact angle of the liquid on the surface. A nearly 180 ° situation is 
depicted in curve (4) of Fig. [H 

For a 100 nm polystyrene (PS) film on a hydrophobized Si wafer with a native oxide 
layer, dewetting starts after heating the sample above the glass transition temperature 
of the polymer. Holes nucleate according to thermal activation or nucleation spots (dust 
particles, inhomogeneities of the substrate or of the polymer film) and grow with time, 
cf. Fig. O The subsequent stages of dewetting are characterized by the formation 
of liquid ridges by coalescence of growing holes and traveling fronts. Very thin films 
in the range of several nanometers may become unstable and can dewet according to 
thermally induced capillary modes, that are amplified by forces contributing to the 
effective interface potential. This phenomenon is characterized by the occurrence of 
a preferred wavelength and is called spinodal dewetting. To study thin film fiow with 
regard to the infiuence of slippage, especially nucleated holes enable an easy experimental 
access for temporal and spatial observation. 

While fronts retract from the substrate and holes grow, a liquid rim is formed at the 
three-phase contact line due to conservation of liquid volume. A common phenomenon 




Figure 2. Dewetting of a 80 nm PS(65 kg/mol) film at T = f35°C from a 
hydrophobized Si substrate captured by optical microscopy (adapted from jlOj). 



is the formation of liquid bulges and so-called "fingers" due to the fact that the liquid 
rim becomes unstable, similar to the instability of a cylindrical hquid jet that beads up 
into droplets (e.g. in case of a water tap). This so-called Ray leigh- Plateau instability 
is based on the fact that certain modes of fluctuations become amplified and surface 
corrugations of a characteristic wavelength become visible. If two holes coalesce, a 
common ridge builds up that in the end decays into single droplets due to the same 
mechanism. The final stage is given by an equilibrium configuration of liquid droplets 
arranged on the substrate exhibiting a static contact angle. Actually, the final state 
would be one single droplet, since the Laplace pressure in droplets of different size is 
varying. Yet then, a substantial material transport must take place, either over the gas 
phase or via an equilibrium film between the droplets. This phenomenon is also called 
Ostwald ripening. In polymeric liquids, these transport pathways are usually extremely 
slow so that a network-like pattern of liquid droplets is already termed "final stage". 
Besides the dynamics and stability of thin liquid films driven by intermolecular forces, 
a recent article by Craster and Matar reviews further aspects such as e.g thermally or 
surfactant driven flows [llj. 

In the first part of this topical review, basic concepts from hydrodynamic theories in 
the bulk situation to corresponding models with regard to confinement are introduced. 
Polymers are regarded as ideal model liquids due to their low vapor pressure, the 
available chemical pureness and furthermore the fact that their viscosity can be 
controlled in a very reliable manner. By setting the viscosity via temperature, the 
experimental conditions can be tuned so that dewetting dynamics can be easily captured. 
Mass conservation can safely be assumed, which consequently simplifies the theoretical 
description. The dewetting dynamics governed by the driving forces and the mechanisms 
of energy dissipation will be discussed also with regard to the shape of hquid ridges. The 
second part summarizes experimental studies concerning the dynamics of two different 
dewetting geometries: the straight front geometry and the growth of holes. Especially 
the influence of parameters such as dewetting temperature, viscosity and molecular 
weight of the polymer will be discussed in detail. Moreover, we focus on various scenarios 
at the solid/liquid interface on the molecular level. Simulations such as for example 
molecular-dynamic (MD) studies help to obtain more information and are supportive 
to gain insight into the molecular mechanism of slippage. 
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2. Basic theoretical concepts 

In this section we aim to describe the main concepts of fluid dynamics especially in 
a conflned geometry. Depending on the type of liquid, viscous or even viscoelastic 
effects have to be considered and deviations from Newtonian behavior might become 
non- negligible. Since we concentrate in this article on polymer melts, viscosity and 
viscoelasticity can be varied by chain length and branching of the polymer. An 
important aspect of a moving liquid is the velocity proflle. 

2.1. Polymer properties 

For a comprehensive understanding of slippage, some important polymer properties 
like glass transition temperature, viscosity and viscoelasticity have to be taken into 
account. Especially in geometries like in a liquid fllm, confinement effects are a concern. 
A detailed description can be found in textbooks [121 US] • 

2.1.1. Polymer physics Polymers are synthesized by polymerization of monomers of 
molar mass M^ono- To characterize the polydispersity of polymer chains in a solution 
or in a melt, the polydispersity index M^/M„ is calculated as the ratio of mean values 
given by the weight and number averaged M„ molecular weight. Polymers are 
able to change their conformations. The radius of gyration characterizes the spatial 
dimension of the polymer and is given as the mean square displacement between 
monomers and the polymer's center of mass. In an isotropic configuration, the shape of 
the polymer chain can be approximated as a spherical entity. Polystyrene, abbreviated 
PS, which is commonly used as a melt in dewetting experiments, is a linear homopolymer 
with Mmono = 104g/mol. Besides other properties concerning the micro-structure of 
polymer chains such as the tacticity and the architecture (linear, branched, ring-shaped), 
physical properties are of special interest. If a melt is heated above its glass transition 
temperature Tg, a phase transition from the glassy phase to the liquid phase occurs 
and the polymer becomes liquid. Randomly structured macromolecules such as atactic 
polymers avoid the formation of semi-crystalline domains below Tg and exhibit a pure 
amorphous phase. The glass transition of a bulk polymer or of a polymeric thin film 
can be observed e.g. via probing its linear expansion coefficient. Although, the glass 
transition is based on a kinematic effect and does not occur due to a rearrangement 
process of polymer chains. Therefore, the glass transition usually takes place on a 
specific temperature range and not at a exactly allocable temperature. According to 
the increased mobility of shorter chains, their glass transition temperature is decreased 
significantly. For PS with a sufficiently large chain-length, Tg=100°C. 

The viscosity ?7 of a polymer melt measures the inner friction of polymer chains 
and governs the time scale of fiow processes. Due to the fact that the mobility of chains 
increases while the temperature increases, the viscosity decreases. Internal stresses relax 
and dynamical processes proceed faster. The most common description of the functional 



Slip effects in polymer thin films 



6 



dependency of the viscosity 77 of a polymer on temperature T was developed by Williams, 
Landel and Ferry: 



In this so-called WLF-equation, rig denotes the viscosity at T^,, B an empirically obtained 
constant, Too the so-called Vogel temperature and fg the free liquid volume fraction. 

Besides the above discussed impact of temperature on the viscosity of a polymer 
melt, also the molecular weight strongly influences rj. While My^ increases, the 
chain mobility decreases and therefore the relaxation times and thus r] increase. For 
sufficiently small M^, the Rouse model predicts a linear increase of rj for increasing M^. 
At this point, another characteristic number of polymer physics has to be introduced: 
the critical chain length for entanglements Mc- Above Mc, the polymer chains form chain 
entanglements exhibiting a specific mean strand length called entanglement length M^. 
For PS, Mc=35kg/mol and Me=17kg/mol is found |12]. According to the reptation 
model of de Gennes [H], the viscosity increases stronger in the presence of chain 
entanglements and an algebraic behavior of rj oc is expected. Empirically, the linear 
regime below Mc is well reproduced, above Mc an exponent of 3.4 is found experimentally 
for different polymers. For a detailed description of chemical and physical properties of 
polymers as well as for the Rouse or the reptation theory we refer to the book Polymer 
physics by Rubinstein and Colby [T2] . 

2.1.2. Properties in confined geometries In contrast to the bulk situation, where 
volume properties of the polymers are measured, liquids in confined geometries such as 
a thin film often show deviations from the behavior in the volume due to additional 
interface effects. One of these properties is the aforementioned glass transition 
temperature. It has been shown in numerous studies, that Tg changes with film thickness 
h. On the one hand, in case of free-standing or supported films exhibiting no or 
repulsive interactions with the substrate, Tg{h) decreases with decreasing film thickness 
[T5| [T6| \rf\ [T8] . On the other hand, Keddie and Jones have shown that an increase of 
the glass transition temperature with decreasing film thickness is possible for attractive 
interactions between substrate and polymer film [19]. The infiuence of the interfacial 
energy on the deviation of Tg from its bulk value has been studied and quantified by 
Fryer and co-workers for different polymers [20]. In case of PS on a solid support, a 
significant change of Tg is found for films thinner than 100 nm (see Fig. [3]) [21]. For 
PS(2k) below 10 nm for instance, the glass transition temperature and the viscosity of 
the polymer film are affected such that these films are liquid at room temperature and 
may dewet spinodally. 

Several attempts have been made to explain the change of Tg according to the 
film thickness. Besides interface-related effects such as reorientation of polymer chains 
or accumulation of chain-ends at the interface, finite-size effects have been proposed 
to be responsible. Herminghaus et al discussed the strain relaxation behavior of thin 



rj = rjg exp 



B{Tg-T) 
fgiT-T^y 



(2.1) 
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Figure 3. Glass transition temperature Tg of polystyrene films of 2kg/mol against 
film thickness (adapted from [21). 

viscoelastic polymer films with regard to surface melting and the shift of the glass 
transition temperature [21]. Kawana and Jones studied the thermal expansivity of 
thin supported polymer films using ellipsometry and attributed their results concerning 
Tg to a liquid-hke surface layer [22], a result that was also found by other authors 
[23| [21]. Besides confinement effects on Tg, further interface-related phenomena have 
been studied: Si et al have shown that polymers in thin films are less entangled than 
bulk polymers and that the effective entanglement molecular weight Mg is significantly 
larger than the bulk value [25] . 

2.1.3. Viscosity and viscoelasticity One of the major characteristics of a liquid in 
general or a polymer in particular is its viscosity rj. Applying shear stress a to a liquid, 
it usually reacts with a strain 7. If stress and strain rate 7 are proportional, the fluid is 
called Newtonian. The constant of proportionality is identified as the viscosity rj of the 
liquid. 

a = r/7. (2.2) 

Liquids such as long-chained polymers show a shear rate dependent viscosity 77(7) 
due to the fact that the liquid molecules are entangled. If the viscosity increases while 
shearing the liquid, we call this behavior shear thickening, whereas in case of lowered 
viscosity so-called shear thinning is responsible. In contrast to the elastic deformation 
of a solid, a deformation of a viscoelastic liquid might induce an additional flow and can 
relax on a specific time scale r. On short time scales (t < r), the liquid behaves in an 
elastic, on long time scales (t > r) in a viscous manner. Thereby, strain 7 is connected 
to stress a via the elastic modulus G of the hquid: 



cr = G'y. 



(2.3) 
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Figure 4. Maxwell model represented by a dashpot and a spring in a serial connection. 

To cover the stress relaxation dynamics of a polymer film, several modelling 
attempts have been proposed. Mostly, so-called Maxwell or Jeffreys models are applied. 
The simplest model is the Maxwell model (see Fig. Hj), which assumes a serial connection 
of a perfectly elastic element (represented by a spring) and a perfectly viscous one 
(represented by a dashpot). Consequently, the total shear strain 7 is given by the sum 
of the corresponding shear strains 7e and 7^ of both mechanisms: 

7 = 7e + Iv (2.4) 

With (El and (E^D we get 

cr = Gule = VmIv (2.5) 

since both react on the same shear stress a. Thereby, the ratio of viscosity of the viscous 
element to the elastic modulus of the elastic one can be identified with a specific time 
scale, the relaxation time tm'- 

Tm = (2.6) 

The relaxation of stress after a step strain 7 leads to a time-dependent stress function 
a{t) for a viscoelastic liquid. Due to the fact that the total strain 7 is constant, a first 
order differential equation for the time-dependent strain 7«(t) is obtained: 

tmIv = 1 - lv{t)- (2.7) 

Solving this differential equation using the initial condition 7^(t = 0) = gives a simple 
exponential decay of a{t) on the time scale of the stress relaxation time tm in the 
Maxwell model: 

7e(t) = 7exp {-t/TM), cr(t) = Guleit) = GA/7exp {-t/rM)- (2.8) 

A situation of special interest is the linear response region: For sufficiently small values 
of 7, the stress/strain-relation (12.31) is valid and the stress relaxation modulus G(t) is 
independent of the strain 7. In this regime, a linear superposition of stresses resulting 
from an infinite number of strain steps can be used to model a steady simple shear fiow 
of a viscoelastic liquid. For larger applied shear rates, linear response and the linear 
superposition fails. The viscosity is still defined as the ratio of stress and strain rate, 
but it has to be regarded as an apparent viscosity which differs from the above described 
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"zero shear rate" viscosity. Polymers with shear-thinning or shear-thickening properties 
can be described by the function 

a oc 7", (2.9) 

where the exponent n can be extracted from experimental data. These type of fluids are 
also called "power law fluids". Moreover, also other nontrivial stress-strain relations can 
be considered or alternatively non-linear extensions can be applied to the linear Maxwell 
models. In case of the linear Jeffreys model, the stress tensor aij relaxes according to 
the following constitutive relaxation equation: 

{l + \idt)a = 7]{l + \2dt)^, (2.10) 

where the strain rate is given by the gradient of the velocity field = djUi + diUj. 
Hereby, Ai governs the relaxation of stress, whereas A2 (A2 < Ai) describes the 
relaxation of the strain rate, respectively. This model accounts for the viscous and 
the elastic properties of a fluid and was used by Blossey, Rauscher, Wagner and Miinch 
as basis for the development of a thin-film equation that incorporates viscoelastic effects 
p6t [27l [28] . For a more elaborate description of non-Newtonian flows we refer e.g. to 
the correspondent work of Te Nijenhuis et al [29]. 



2.1.4. Reynolds and Weissenberg number The flow of a liquid can be characterized by 
specific numbers. One of these numbers is the so-called Reynolds number Re, which 
describes the ratio of inertia effects to viscous flow contributions. In case of thin liquid 
films. Re can be written as 

i?e=^, (2.11) 
V 

where p denotes the density of the liquid, u describes the flow velocity and h stands 
for the film thickness [30]. For thin dewetting polymer films, the Reynolds number is 
very small, i.e. Re <C 1, and a low-i?e lubrication theory can be applied. To quantify 
and to judge the occurrence of viscoelastic effects versus pure viscous flow, the so-called 
Weissenberg number Wi has been introduced as 

Wi = Tj. (2.12) 

Thereby, r denotes the relaxation time and 7 the strain rate as introduced in the previous 
section. If Wi <^ 1, an impact of viscoelasticity on flow dynamics can be neglected and 
viscous flow dominates. 



2.2. Navier-Stokes equations 

The Navier-Stokes equations for a Newtonian liquid mark the starting point for the 
discussion of fluid dynamics in confined geometries. According to conservation of mass, 
the equation of continuity can be formulated as 



dtp + V- (pu) = 0, 



(2.13) 



Slip effects in polymer thin films 



10 



where u = {ux,Uy,Uz) is the velocity field of the fiuid. For an incompressible liquid, 
which implies a temporally and spatially constant liquid density p, f l2.13p can be 
simplified to 

Vu = 0. (2.14) 

With the conservation of momentum, the Navier-Stokes equations for an incompressible 
liquid can be written as 

p(9t + u ■ V)u = -Vj9 + T^Au + f, (2.15) 

with the pressure gradient Vp and the volume force f of external fields acting as driving 
forces for the liquid fiow. We already stated that for small Reynolds numbers, i.e. 
i?e <^ 1, the terms on the left hand side of (12.151) can be neglected as compared to terms 
describing the pressure gradient, external volume forces and viscous fiow. By that, we 
can simplify (I2.15P to the so-called Stokes equation 

= - Vp + r/Au + f. (2.16) 

In section 12.51 we will demonstrate how these basic laws of bulk fiuid dynamics can be 
applied to the fiow geometry of a thin film supported by a solid substrate. 



2.3. Free interface boundary condition 

At the free interface of a supported hquid film, i.e. at the liquid/gas or usually the 
liquid/air interface, no shear forces can be transferred to the gas phase due to the 
negligible viscosity of the gas. In general, the stress tensor a*j is given by the stress 
tension aij, see (12. 2p . and the pressure p: 

a*j = (Tij + p6ij = r]{djUi + diUj) + p6ij. (2.17) 

The tangential t and normal n (perpendicular to the interface) components of the stress 
tensor are: 

(cT* ■ n) ■ t = ((T* ■ n) ■ n = 7,„k, (2.18) 

where k denotes the mean curvature and 7/^, the interfacial tension (i.e. the surface 
tension of the liquid) of the liquid/vapor interface. If the liquid is at rest, i.e. the 
stationary case u = 0, the latter boundary condition gives the equation for the Laplace 
pressure pl- 

11 , , 

PL = 7lvK = -flv{-^ + -^)- (2.19) 

Ri and R2 are the principal radii of curvature of the free liquid/gas boundary; the 
appropriate signs of the radii are chosen according to the condition that convex 
boundaries give positive signs. Such convex liquid/gas boundaries lead to an additional 
pressure within the liquid due to its surface tension. In the next section, the solid/liquid 
boundary condition will be discussed, which yields a treatment of slip effects. 
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Figure 5. Different velocity profiles in the vicinity of the solid/liquid interface and 
illustration of the slip (extrapolation) length b. The situation of so-called "apparent" 
slip is illustrated on the right: According to a thin liquid layer of thickness zq 
that obtains a significantly reduced viscosity, the slip velocity Ux\z=o is zero, but a 
substantial slip length is measured. 



2.4- Slip/no-slip boundary condition 

2.4-1 ■ Navier slip boundary condition In contrast to fluid dynamics in a bulk volume, 
where the assumption that the tangential velocity U|| at the solid/liquid interface 
vanishes [no-slip boundary condition), confined geometries require a more detailed 
investigation as slippage becomes important. In 1823, Navier [21] introduced a linear 
boundary condition: The tangential velocity U|| is proportional to the normal component 
of the strain rate tensor; the constant of proportionality is described as the so-called 
slip length h: 



u 



6n ■ 7 



(2.20) 



In case of simple shear flow in a;-direction, the deflnition of the slip length can be 
alternatively written as 



2=0 



a 



(2.21) 



where ^ = a/ux denotes the friction coefficient at the solid/liquid interface. The xy- 
plane thereby represents the substrate surface. According to these definitions, the slip 
length can be illustrated as the extrapolation length of the velocity profile " inside" the 
substrate, cf. Fig. [5l Moreover, both limiting cases are included within this description: 
For 6 = 0, we obtain the no-slip situation, whereas b = oo characterizes a full-slip 
situation. The latter case corresponds to "plug-fiow", where the liquid behaves like a 
solid that slips over the support. 



2.4-2. How to measure the slip length? In recent years, numerous experimental studies 
were published using diverse methods to probe the slip length at the boundary of 
different simple or complex liquids and solid supports. For details concerning these 
experimental methods we refer to the review articles from Lauga et al [32], Neto et 
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al [33] and Bocquet and Barrat [31] (and references therein). To probe the boundary 
condition, scientists performed either drainage experiments or direct measurements of 
the local velocity profile using e.g. tracer particles. 

In case of drainage experiments, the liquid is squeezed between two objects, e.g. a 
flat surface and a colloidal probe at the tip of an AFM cantilever, and the corresponding 
force for dragging the probe is measured {colloidal probe AFM). Alternatively, in an 
surface force apparatus ( SFA ), two cylinders arranged perpendicular to each other are 
brought in closer contact and force/distance measurements are performed to infer the 
slip length. 

The use of tracer particles as a probe of the local flow profile might bring some 
disadvantages. The chemistry of these particles is usually different from the liquid 
molecules and their influence on the results might not be negligible. A similar method 
is called fluorescence recovery after photo bleaching. Thereby, a distinct part of a 
fluorescent liquid is bleached by a laser pulse and the flow of non-bleached liquid into 
that part is measured. The disadvantage of this method is that diffusion might be a 
further parameter that is hard to control. Recently, Joly et al showed that also thermal 
motion of confined colloidal tracers in the vicinity of the solid/liquid interface can be 
used as a probe of slippage without relying on external driving forces [35] . 

2.4-3. Which parameters influence slippage ? Of course, many interesting aspects in the 
field of micro- and nanofiuidics are related to intrinsic parameters that govern slippage 
of liquid molecules at the sohd/liquid interface. For simple liquids on smooth surfaces, 
the contact angle is one of the main parameters influencing slippage [361 EH EH [39] . 
This originates from the effect of molecular interactions between liquid molecules 
and the solid surface: If the molecular attraction of liquid molecules and surface 
decreases (and thereby the contact angle increases), slippage is enlarged. Further 
studies aim to quantify the impact of roughness [371 SOI SIl 112] or topographic structure 
[43l l44l l45t l46t 147] of the surface on slippage. For different roughness length scales, 
a suppression (see e.g. [371 SO]) or an amplification (see e.g. ^3]) of slippage can be 
observed. Moreover, the shape of molecular liquids itself has been experimentally shown 
to impact the boundary condition. Schmatko et al found significantly larger slip lengths 
for elongated linear compared to branched molecules [I8j|. This might be associated 
with molecular ordering effects [19] and the formation of layers of the fluid in case of the 
capability of these liquids to align in the vicinity of the interface [50]. Cho et al identified 
the dipole moment of Newtonian liquids at hydrophobic surfaces as a crucial parameter 
for shp [51]. De Gennes proposed a thin gas layer at the interface of sohd surface and 
liquid as a possible source of large slip lengths [52]. Recently, MD studies for water on 
hydrophobic surfaces by Huang et al revealed a dependence of slippage on the amount of 
water depletion at the surface and a strong increase of slip with increasing contact angle 
|53j . Such depletion layers for water in the vicinity of smooth hydrophobic surfaces 
have also been experimentally observed using scattering techniques [HH [551 [HEl [57] . 
Contamination by nanoscale air bubbles (so-called nanobubbles) and its influence on 
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Figure 6. Illustration of the nomenclature of the thin film length scales (x and y are 
parallel to the substrate) and the velocity contribution U|| ~ {ux,Uy). 

slippage has been controversially discussed in literature (see e.g. [58|, [59l [601 EI])- In 
case of more complex liquids such as polymer melts further concepts come into play. 
They will be illustrated in section 13. 4[ 

2. 5. Thin-film equation for Newtonian liquids 

2.5.1. Derivation Confining the flow of a liquid to the geometry to the one of a 
thin film, we can assume that the velocity contribution perpendicular to the substrate 
is much smaller than the parallel one. Furthermore, the lateral length scale of film 
thickness variations is much smaller than the film thickness itself. On the basis of these 
assumptions, Oron et al [62] developed a thin-film equation from the rather complex 
equations of motion, (12.131) and (I2.15p . In case of film thicknesses smaller than the 
capillary length = \/jtJ~pg, (which is typically in the order of magnitude of 1 mm) 
fl2.16p can be written as 

= -V(p + 0'(/i)) +r7Au (2.22) 

Additional external fields such as gravitation can be neglected, but a secondary 
contribution 4>'{h), the disjoining pressure, has been added to the capillary pressure 
p. The disjoining pressure originates from molecular interactions of the fluid molecules 
with the substrate. The effective interface potential 0(/i) summarizes the inter-molecular 
interactions and describes the energy that is required to bring two interfaces from infinity 
to the finite distance h. As already discussed in the introductory part, the stability of a 
thin liquid film is also governed by (j){h). For a further description of thin film stability, 
we refer to [8] and the references therein. 

The derivation of a thin-film equation for Newtonian liquids starts with the 
kinematic condition 

d^h = -V|| [ u\\dz, (2.23) 
Jo 

i.e. the coupling of the time derivative of h{x, y, t) to the flow field, where the index 
II in general denotes the components parallel to the substrate (V|| = {dx,dy) and 
U|| = {ux,Uy)) as illustrated in Fig. [61 
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For thin liquid films, film thickness variations on lateral scale L are much larger 
than the length scale of the film thickness H. Introducing the parameter e = H/L ^ 1 
yields the so-called lubrication approximation and is used in the following to re-scale 
the variables to dimensionless values. In a first approximation, linearized equations are 
obtained while neglecting all terms of the order O(e^). For reasons of simplicity and due 
to translational invariance in the surface plane, a one-dimensional geometry is used: 

d^{p + (t)') = dlu^, 9,(p + 0')=O, d,u, + d,u, = Q. (2.24) 

While the substrate is supposed to be impenetrable for the liquid, i.e. = for 2; = 0, 
friction at the interface implies a velocity gradient d^u^ = u^/b for z = 0. Moreover, 
the tangential and normal boundary condition at the free interface, i.e. z = h{x), can 
be simplified in the following manner: 

d,u^ = 0, p + dlh = 0. (2.25) 

From ( 12.24P and the boundary conditions, the velocity profile Ux{z) can be obtained. 
Using the kinematic condition f l2.23p . the equation of motion for thin films in three 
dimensions is derived: 

d,h = -V[m{h)V{^i,Ah - <f)'{h))], (2.26) 

where m{h) denotes the mobility given by 

m(h) = —{h'' + 3bh^). (2.27) 

2.5.2. Lubrication models including slippage As discussed in the previous section, the 
derivation of the thin-film equation is based on the so-called lubrication approximation 
and the re-scaling of relevant values in e. As a consequence, the slip length b is supposed 
to obtain values smaller than the film thickness /i, i.e. b <^ h. To extend this so-called 
weak-slip situation with regard to larger slip b^ Miinch et al [30] and Kargupta et 
al [63] developed independently so-called strong-slip models. Thereby, the slip length is 
defined as 6 = /?/e^. The corresponding equation of motion together with the kinematic 
condition in one dimension for a Newtonian thin liquid film read as: 

u = —dx{2rihdxu) + —dx{'-fiydlh-(j)'{h)) 

V V (2.28) 

dth = -d^iyhu). 

In fact, a family of lubrications models, cf. Tab. [H accounting for different slip 
situations have been derived. In the limit 6 — 0, i.e. the no-slip situation, the mobility 
is given by m{h) = h^/3r]. If the slip length is in the range of the film thickness b h, 
the mobility in the corresponding intermediate- slip model is m{h) = bh'^/rj. Recently, 
Fetzer et al [64] derived a more generalized model based on the full Stokes equations, 
developed up to third order of a Taylor expansion. The authors were able to show that 
this model is in good agreement with numerical simulations of the full hydrodynamic 
equations and is not restricted to a certain slip regime as the aforementioned lubrication 
models. 
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Table 1. Summary of lubrication models for Newtonian flow and different slip 
situations. 



model 


validity 


equation 


limiting cases ref. 


weak-slip 


b^h 


(12.261). (12.2711 


b^O (no-slip) 

6 — > oo (intermediate-slip) 


strong-slip 


b:> h 


(|2.28D 


/? -> (intermediate-slip) _ _ 
/?->oo ("free" -slip) 



2.5.3. Lubrication models including viscoelasticity In the meanwhile, the derivation of 
a thin film equation for the weak-slip case including linear viscoelastic effects of Jeffreys 
type (such as described by equation fl2.10p in section [2. 1.3p has been achieved (see [2^]). 
To cover relaxation dynamics of the stress tensor cr, an additional term V ■ cr on the 
right hand side of (I2.22p has to be included to the aforementioned model for Newtonian 
liquids. Furthermore, the treatment of linear viscoelastic effects was also achieved for 
the strong-slip situation by Blossey et al [27]. To summarize these extensions, the 
essential result is the fact that linear viscoelastic effects are absent in the weak-slip case 
and the Newtonian thin-film model is still valid. The strong-slip situation, however, 
is more complicated. Shppage and viscoelasticity are combined and strongly affect the 
corresponding equations. In the meanwhile, the authors were able to fully incorporate 
the non-linearities of the co-rotational Jeffreys model for viscoelastic relaxation into 
their thin-film model [28] . 

The extensions of the aforementioned thin-film models for different slip conditions 
with or without the presence of viscoelastic relaxation (Newtonian and non-Newtonian 
models) affect on the one hand the rupture conditions, but also on the other hand 
the shape of a liquid ridge. These two phenomena will be discussed in the next two 
subsections. A elaborate description of these theoretical aspects can be found in a 
recent review article by Blossey [65j . 

2.5.4. Application I - Spinodal dewetting One of the main applications of the 
theoretical thin-film models is the dewetting of thin polymer films. As introduced in 
section [1] and illustrated by Fig. [1], the stability of a thin liquid film is governed by the 
effective interface potential. Basically, long-range attractive van der Waals forces add 
to short-range repulsive forces. Due to the planar geometry of two interfaces of distance 
d, the van der Waals contribution to the potential is (f>{d)vdw oc —A/d"^, where A is the 
Hamaker constant. For the description of the explicit calculation of Hamaker constants 
from the dielectric functions of the involved materials we refer to [Bj and to the book by 
Israelachvili [66]. Experimental systems often exhibit multi-layer situations, cf. Fig. [71 
A hydrophobic film and/or an oxide layer of distinct thicknesses di exhibiting Hamaker 
constants Ai require a superposition of contributions to the potential: 
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Figure 7. Two examples for polystyrene films prepared on multi-layer substrates. 
Left: Silicon wafer with native oxide layer covered with a hydrophobic layer (e.g. a 
self- assembled monolayer (SAM)). Right: Silicon wafer with an increased (compared 
to a native silicon oxide) oxide layer thickness. 




Figure 8. Effective interface potential (pih), gained from experimental data, plotted 
against the thickness /i of a thin PS film prepared on silicon wafers with different oxide 
layer thicknesses di (adapted from |/7j). The cross-hatched rectangle marks the error 
for heq and the depth of the global minimum. 



Consequently, the shape of the effective interface potential (and thereby also the thin 
film stability) is governed by the set of Hamaker constants Ai and film thicknesses d^. 
E.g. for a thin PS film of thickness /i on a Si substrate with a native oxide layer of 
2.4 nm, (j){h) shows a global minimum at an equilibrium film thickness h^q (c.f. Fig. [8]). 
Moreover, a local maximum ai h > hf.q is found. 

If the PS film is sufficiently thin {(f)" (h) < 0) , it may become unstable due to 
the amplification of thermally induced capillary waves. To track the evolution of small 
fiuctuations of the film thickness h, i.e. /(x, t) = /i(x, t) — h with /(x, t) <^ h, a. 
Fourier transform of the linearized thin film equation (I2.26P has to be performed. The 
amplitudes of the capillary waves grow exponentially in time. Their growth rate a can 
be calculated as a function of the wavenumber q and depends on the sign of the local 
curvature of the interface potential. If the second derivative of the effective interface 
potential (p" at the film thickness h is positive, a is negative for all q and the amplitudes 
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Figure 9. In situ (at elevated temperature T = 53 °C) atomic force microscopy (AFM) 
images (corresponding annealing times given in pictures) of a spinodally dewetting 
3.9(2) nm PS(2kg/mol) film on a Si wafer with a thick (191 nm) oxide layer (adapted 
from [67]). 



of the capillary waves are damped. If 0" < 0, the growth rate a is positive for a certain 
range of wavenumbers up to a critical wavenumber qc and capillary waves are amplified. 
The wavenumber that corresponds to the maximum value of a and therefore exhibits 
the fastest amplification, is called preferred wavenumber go and is connected to the 
preferred wavelength Aq = 27r/go- The latter is also called spinodal wavelength and can 
be written as 



The spinodal dewetting process can be monitored e.g. by atomic force microscopy 
(AFM) as shown in Fig. [9l By measuring Aq as a function of film thickness, 4>"{h) 
can be inferred and conclusions can be drawn with regard to the effective interface 
potential (f){h) [7j. For further details concerning stability of thin films we refer to [8] 
and references therein. 

Using the strong-slip model while taking slip into account (12.281) . Rauscher and 
coworkers could theoretically show that slippage is supposed to influence the capillary 
wave spectrum due to a different mobility at the solid/liquid interface [68]. The position 
of the maximum go shifts to smaller wavenumbers and larger wavelengths for increasing 
slip length. As for today, to the best of our knowledge experimental studies concerning 
the impact of slippage on the spinodal wavelength are not available. 
The description of the influence of thermal noise on the temporal and spatial dynamics 
of spinodally dewetting thin polymer films has been recently achieved by Fetzer and 
coworkers [69]. A stochastic Navier-Stokes equation with an additional random stress 
fluctuation tensor that accounts for thermal molecular motion is utilized to model the 
flow while assuming a no-slip boundary condition at the solid/liquid interface. The 
stochastic model matches the experimentally observed spectrum of capillary waves and 
thermal fluctuations cause the coarsening of typical length scales. 




(2.30) 
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Figure 10. Left: AFM image (scan size 10 /xm) of a liquid rim formed during hole 
growth in a PS film dewetting from a hydrophobic substrate (Si wafer covered with 
a 21 nm hydrophobic Teflon® coating (AF1600)). Right: Different rim morphologies 
for PS on AF 1600 (AFM cross-sections): Profile exhibiting a trough (depicted in the 
inset) for PS(35.6 kg/mol) at 120 °C and a monotonically decaying rim shape for 
PS(125kg/mol) at 130 °C. 

2. 5. 5. Application II - Shape of a dewetting rim Besides the implications on spinodal 
dewetting, the one-dimensional thin-film model has been successfully applied to the 
shape of the rim along the perimeter of e.g. nucleated holes. Experimentally, researchers 
have studied and observed different types of rim profiles [TOl [TH [72]. As shown in Fig. 
[TOl profiles either decay monotonically into the undisturbed film or they show a more 
symmetrical profile exhibiting a trough (termed "oscillatory shape"). If the depth of 
this through is in the range of the film thickness, a ring of so-called satellite holes can 
be generated [73l [71] . 

The shape of a dewetting rim can be understood by the aforementioned thin film 
theory for Newtonian liquids: Introducing a small perturbation 5h{x,t) -C h of the film 
thickness h{x,t) and small velocities u{x,t) leads to linearized thin-film equations that 
describe the temporal and spatial evolution of 6h. Thereby, the disjoining pressure (p'{h) 
can be neglected due to the fact that films thicker than 10 nm are considered. To obtain 
stationary solutions of the linearized equations, a frame that is co-moving with the rim 
C{x,t) = X — s{t) is introduced. Thereby, s{t) denotes the position of the three-phase 
contact line; s stands for the dewetting velocity V as described in the next sections. 
Fetzer et al used a normal modes ansatz 6h{Q = 6hoexpk( and u{() = uoexpk( in 
the linear stability analysis, which leads to a characteristic polynomial of third order. 
Depending on the ratio of slip length to film thickness b/h and on the capillary number 
Ca = ris/jiy, the parameter k obtains complex or real solutions. Fetzer and coworkers 
successfully identified the morphological transition from oscillatory to monotonically 
decaying rims and were able to extract slip lengths and capillary numbers from diverse 
experiments on dewetting surfaces [751 [13 [HH [77] . 
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Figure 11. Temporal series of optical micrographs showing the growth of a hole in 
a 120nm PS(13.7kg/mol) film at T = 120 °C prepared on a Si wafer covered with 
a 21 nm hydrophobic Teflon® coating (AF1600). In the last stage of hole growth 
(right image) perturbations of the three-phase contact line (according to the liquid 
rim instability) become visible. 



3. Flow dynamics of thin polymer films - Experimental studies and 
theoretical models 

One of the main aspects of experimental studies concerning the flow dynamics of thin 
polymer fllms is to obtain a comprehensive view on the molecular mechanisms of slippage 
and on the responsible parameters. Although these insights are rather indirect, several 
models have been proposed to explain diverse experimental results. In this section, we 
will focus on these studies, with special regard to the proposed mechanisms of slippage 
at the solid/polymer interface. In general, we have to distinguish two different dewetting 
geometries: growth of holes (cf. Fig. [TTj) and receding straight fronts. On the one hand 
driving forces and on the other hand dissipation mechanisms have to be considered. 

3.1. Dewetting dynamics - Driving forces 

According to the description of the effective interface potential in section 12.5.41 a global 
minimum of occurs at h^.^ in case of unstable or metastable fllms. This means that 
the fllm will thin until a thickness of heq is reached. In other words, a thin wetting 
layer remains on top of the substrate, if heq > and if heq has a size that is not below 
the size of the molecules 0. After dewetting has taken place (the dynamics of which 
is not covered by </)(/i) but depends on viscosity, viscoelasticity, contact angle and the 
solid/liquid boundary condition), single droplets remain on top of the wetting layer. 
The droplets - in equilibrium - exhibit the Young's contact angle 6y- Parallel to the 
Young's equation that characterizes the contact angle via the involved surface tensions, 
another characteristic parameter can be deflned to specify the wettability of a surface 
by a liquid. This is the so-called spreading parameter S given by 

S = Jsv - ilsl + llv + (p{heq)), (3.1) 

I Note that the concept of the effective interface potential is a continuum approach. It may fail if 
molecular sized phenomena are to be captured, such as extremely thin films close to the size of the 
molecules. 
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where denotes the interface tension of the sohd (s), hquid (/) and the vapor {v) 
phase. S describes the energetic difference per unit area between a dry surface and a 
surface that is covered by a hquid layer. Consequently, if < the system can gain 
energy by reducing the contact area with the substrate. Thereby, the liquid forms a 
dynamic contact angle 9d = 9y/\/^ at the three-phase contact line [HT] . 

The link between the contact angle in equilibrium and the depth of the minimum of 
the effective interface potential becomes obvious, if the equation for the Young's contact 
angle is inserted in fl3.1l) : 

(t){heq) = 'Jlvicosdy - I) = S (3.2) 

Driving forces for wetting and dewetting are capillary forces (see [ZH])- Thereby, 
the capillary force per unit length of the contact line is in general given by 
Fc 

-J = (cos 6'y - cosOd). (3.3) 

As long as 6^ < Oy remains valid, the capillary force is negative and the three-phase 
contact line will recede; the liquid film dewets. Dewetting ends as soon as droplets 
exhibit their Young's contact angle Oy on the surface. During dewetting, a rim which 
is formed from accumulated material and retracts from the substrate while further 
growing. De Gennes developed a theoretical description of the resulting driving force 
for dewetting [78] : 

On the "dry" side of the rim, which denotes the side where the three-phase contact 
line is located, a negative capillary force pulls at the contact line: 

^ = S + -ii,{l- cosed). (3.4) 

On the other side, where the rim decays into the liquid film (the "wet" side of the rim), 
a positive capillary force occurs: 
Fc 

-J = '^i^{l-cos9d). (3.5) 

We have to remark that in experimental systems, the Laplace pressure of course avoids 
edges and therefore forces the dry side of the rim to decay smoothly into the film. 
Summing up both forces per unit length gives the effective driving force for retracting 
contact lines in dewetting systems: 

y = 1^1 = 7«.(1 - COS^^y) ~ ^li^Oy. (3.6) 

Thereby, it is assumed that the system is in a quasi-stationary state, which means that 
changes in the shape of the rim occur much slower that the dewetting velocity. All in 
all, this result implies that the driving force for dewetting is given by the absolute value 
of the spreading parameter, which solely depends on the surface tension 7^ of the liquid 
(a property which is purely inherent to the liquid) and the Young's contact angle of 
the liquid on the surface. The same conclusion is drawn if considering energies instead 
of driving forces. The depth of the minimum of the effective interface potential 0(/i) 
is equal to S and gives the energy per unit area that is set free during dewetting [79] . 
Thereby, S also stands for the force per unit length of the contact line (see (13. 6p ). 
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3.2. Dewetting dynamics - Dissipation mechanisms 

In the previous section, we focused on the driving forces for dewetting. The 
experimentally observed dynamics represents a force balance of driving forces and 
friction forces. The resulting dewetting velocity V is connected to the spreading 
parameter S and the occurring energy dissipation mechanisms (Fj gives the friction 
force per unit length of the contact line) and their corresponding velocity contributions 
Vi via the power balance 

\S\V = Y,F^v^. (3.7) 

i 

In case of dissipation due to viscous friction within the liquid (F^,) and dissipation due 
to friction of liquid molecules at the solid/liquid interface (-F,), i.e. slippage with a 
finite velocity Vg = u\z=o, we end up with the following balance (the index v stands for 
'viscous', s denotes 'slip'): 

\S\V = F,v, + FsVs. (3.8) 



3.2.1. Viscous friction According to the work of Brochard-Wyart et al [80], the no-slip 
boundary condition for fluid flow implies a friction force that is proportional to the liquid 
viscosity and the dewetting velocity v and means that the dissipation is solely due to 
viscous friction within the liquid. The largest strain rates occur in the direct vicinity of 
the three-phase contact line. The dissipation is therefore mainly independent from the 
size of the dewetting rim. However, the flow geometry (given by the dynamic contact 
angle 9d) at the contact line influences viscous dissipation. In case of pure viscous flow 
V = v^, the following expression for the velocity contribution is obtained for viscous 
flow [80]: 

V, = CM^-^. (3.9) 

Cv{Od) denotes the constant of proportionality and displays a measure for the flow field 
in the vicinity of the contact line. According to the description of Redon et al |HI], the 
velocity in case of viscous dissipation strongly depends on 9y and can also be written 

as 

oc —el, (3.10) 

where the constant of proportionality accounts for the flow singularity near the contact 
line. In detail, this constant represents a logarithmic factor that includes i.a. a short- 
distance cutoff and accounts for the fact that viscous dissipation diverges near the 
contact line. Moreover, a small impact of the slip length h on this logarithmic factor 
has been experimentally found [HI] by variation of molecular weight (cf. (13.121) ) by 
Redon ^2] . Consequently, only in the case of a no-slip boundary condition, the constant 
of proportionality in fl3.10p and also Cy{9d) in (13.91) is purely independent of slip. 
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Figure 12. Schematic representation of a rim cross-section illustrating the rim width 
w (usually obtained as the distance between the three-phase contact line and the 
position where the rim height is dropped to 1.1 /i), the distance xq for a given film 
thickness h and hole radius R (or dewetted distance D oi a straight front). (Height 
scale in sketch is exaggerated as compared to the lateral length scale.) 



For growing holes of radius R and dewetting velocity V = dR/dt, integration of 
( 13. 9p gives a linear proportionality of the radius versus dewetting time t: 

Roct (3.11) 

3.2.2. Slippage - Friction at the solid/liquid interface Besides viscous dissipation within 
the rim, de Gennes expected a long-chained polymer film to exhibit an exceptional 
amount of slippage. For entangled polymer melts on a non-adsorbing surface, de Gennes 
stated that the slip length should strongly increase with the chain length of the polymer 
[831: 



iV3 

b = a—, (3.12) 

e 

where a is the size of the monomer, N the polymerization index and the entanglement 
length. In case of dominating slippage, a minor contribution of viscous dissipation is 
expected due to very small stresses within the rim. The corresponding analytical model 
(see [80|[82]) is based on the linear friction force per unit length of the contact line given 
by Fs oc ^Vg. Dissipation occurs along the distance of the solid/liquid interface, where 
liquid molecules are moved over the substrate. In general, this distance is identified as 
the width w of the rim (cf. Fig. [T2ll . 

Consequently, Fg oc w can be assumed. If the slip length is introduced according 
to Navier's model hj b = rj/^ (see f l2.2ip in section [2^ . the slip velocity contribution 
Vs is given by: 

vs = (3.13) 

Let's consider the conservation of mass for a growing hole of radius R, 

7r{R + Xoyh = 27rQ{R + w/2), (3.14) 
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and for a straight dewetting front of dewetted distance D, 

{D + XQ)lh = Ql, (3.15) 

where xo stands for the distance as depicted in Fig. [121 Q denotes the cross sectional 
area of a rim (approximated as a half-circle i.e. Q = 7r(xo/2)^ ^ 7r(ty/2)^), and / 
is the length of a straight front, which cancels out. Based on the assumption of self- 
similar growing rims and the fact that Xo/R <^ 1 and xq/D <^ 1, the width of the 
rim w oc y/h^/R for both geometries. The constant of proportionality, however, which 
we name Cs in the following, depends on the shape of the rim and furthermore on the 
dewetting geometry itself: 



w 



CsVhVR. (3.16) 



Values for Cg can be obtained collecting a temporal series of AFM snapshots of rim 
profiles and fitting fl3.16p to the measured rim width w plotted versus the hole radius 
R. The initial film thickness h can also be measured by AFM. Of course, the validity 
of this simplified model and especially the approximation w ^ xo becomes questionable 
for asymmetric rims and should lead to systematic differences in Cg- We point out that, 
as described in section 12.5.51 the degree of asymmetry of the cross-section of the rim is 
strongly correlated to the capillary number and the ratio of slip length to film thickness. 
Replacing the above derived relation of rim width w and hole radius R in ( 13.13P gives 
the slip velocity contribution in terms of the hole radius R for purely slipping liquids: 

Ms\ b 1 

Vg = ^—=. (3.17) 

3 T] CsVhVR 

Via integration, a characteristic growth law for the radius i? of a dewetting hole with 
time t is found (separation of variables), 

R(xt^'^. (3.18) 

Using the lubrication models described in section 12.51 Miinch compared numerical 
simulations of polymer melts dewetting from hydrophobized substrates to the above 
described dewetting dynamics obtained from scaling arguments (based on energy 
balances) [HI]: In the no-slip situation (mobility m{h) oc h?), an exponent a = 0.913 
is reported instead of 1. The deviation is traced back to the fact that the logarithmic 
factor in the constant of proportionality of equation (I3.10p also depends on the width of 
the rim (which evolves with time t). Numerical simulations of the slip-dominated case 
{m{h) oc h'^) give a = 0.661, which captures a = 2/3 very well. 

§ Note that usually a dynamic contact angle < 90 ° at the three-phase contact line is obtained in 
experiments. The difference between a segment of a circle and the half-circle approximation is marginal 
and contributes to the constant of proportionality Cg- 
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Figure 13. Left: Optical measurement of the hole radius versus time in 130 nm 
PS(13.7kg/mol) films on different substrates. Right: Plot of the dewetting velocity V 
versus l/\fR. The y-axis intercept is identical for all substrates. The difference in the 
slope K indicates substantial differences in slip for dewetting PS (adapted from [55] 
and [87]). 



t - to = - 2§^v^ + ^i^f In (1 + ^)) (3.19) 



3.2.3. Models based on the superposition of dissipation mechanisms Considering the 
models derived before for pure viscous flow and for pure slippage, a combination of 
these models seems to be reasonable for situations, where viscous dissipation as well as 
dissipation at the solid/hquid interface act together. Jacobs et al proposed an additive 
superposition of the corresponding velocity contributions, i.e. V = Vg+v^ [85j, according 
to the fact that both mechanisms counteract the same driving force 15*1. Separation of 
variables leads to the implicit function: 



with 



In that notation, the velocity contributions are given by = \S\/Kv and Vs = 
\S\/{\/RKy). One of the preconditions of this superposition was the successful inclusion 
of both border cases: For b ^ 0, R (x tis obtained, whereas for b oo, R (x t^/'^ follows. 

As a more practicable alternative for fitting this relation between radius R and time 
t to experimental data (typical data is shown in Fig. [T3ll . Fetzer and Jacobs recently 
proposed a more simple way of visualizing slip effects [86]: According to the additive 
superposition, the dewetting velocity can be written in terms of 

V = v^ + ^, ^=^— TT^r- (3.21) 



"^-ckr "'-T^' 1-'^- 



which leads to a linear relationship when plotting the dewetting velocity v versus 1 / y/R 
(c.f. Fig. [T3|) . Then, the y-axis intercept of the straight line can be identified as 
the viscous velocity contribution v^, whereas the the slope K is connected to the slip 
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length as illustrated in fl3.2ip . The validity of this model, which assumes the linear 
superposition of velocity contributions, was checked by plotting the viscous velocity 
contribution versus the reciprocal viscosity r]. An excellent agreement has been 
demonstrated in view of melt viscosity data obtained from independent measurements. 
Moreover, the slope has been used to calculate the slip length b. In the experiments, 
substantial differences in slip lengths for different substrates (cf. Fig. [13]) at identical 
liquid properties could be detected [86] . 

3.2.4- Molecular-kinetic theory and further approaches Besides theoretical models 
based on continuum hydrodynamics (see I3.2.ip . contact line dynamics such as the 
spreading of a liquid droplet on a planar surface has been analyzed according to the 
so-called molecular-kinetic theory (MKT) [HHl EHl EO]. The movement of a contact 
line is described as an activated rate process, where the liquid molecules close to the 
substrate jump from one potential well to another. The resulting contact line friction 
coefficient is proportional i.a. to the viscosity r] of the liquid and to the exponential 
of the reversible work of adhesion 7^(1 + cos^y)/fcBT. As a consequence, increasing 
the temperature reduces the friction coefficient and increases the velocity. In case of 
a squalane droplet spreading on gold surfaces covered with different self-assembled 
monolayers (SAM) of thiols, the corresponding contact line friction coefficient turns 
out to increase linearly with the chain length of the SAM [91]. SAM surfaces have also 
been target of friction experiments using AFM tips. Barrena et al explained discrete 
changes in the frictional behavior with discrete molecular tilts of the chains of the SAM 
[92] . Viscoelastic deformation of the substrate at the contact line due to the normal 
component of the liquid surface tension has also attracted attention with regard to fluid 
dynamics on the nanoscale and energy dissipation. The phenomenon of a reduced 
velocity of contact-line displacements compared to a corresponding non-deformable 
substrate (termed "viscoelastic braking") was extensively studied and described by 
Shanahan and Carre [OH [941 [95l [96l [97]. In this context. Long et al theoretically 
examined the static and dynamic wetting properties of liquids on thin rubber films 
[98] and grafted polymer layers [99]. To the best of our knowledge, there is to-date 
no experimental evidence of deformations and tilting on the molecular scale in case of 
SAM surfaces due to retracting contact-lines of a dewetting polymer film. However, 
dewetting on soft deformable rubber surfaces has recently been shown to be influenced 
by viscoelastic deformations of the substrate [TOOj. In case of no n- volatile liquids such 
as polymers, evaporation into and condensation from the vapor phase can be excluded 
as relevant mechanism of energy dissipation. 

3. 3. Dynamics of growing holes 

Besides very early theoretical and experimental studies of thin film rupture [HI I101[ 
11021 11031 11041 11051 1106] . one of the first studies concerning dewetting experiments had 
been published by Redon et al [81j. Films of alkanes and polydimethylsiloxane (PDMS) 
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have been prepared on different hydrophobized silicon wafers. The authors observed a 
constant dewetting velocity which was inversely proportional to the viscosity and very 
sensitive to the equilibrium contact angle of the liquid on a corresponding surface. In 
1994, Redon et al investigated the dewetting of PDMS films of different thicknesses 
on silanized hydrophobic surfaces [82]. They found that for films larger than 10 /xm a 
constant dewetting velocity is found, whereas R{t) oc t^/^ in films thinner than 1 /xm 
is observed. These results corroborate the models represented by (13. lip and fl3.18p and 
indicate that for sufficiently thick films {h ^ h) viscous dissipation dominates whereas 
slippage becomes relevant in case of thin films {h < h) |82j . 

3.3.1. Stages of Dewetting While studying the growth of holes, Brochard-Wyart et 
al proposed a set of subsequent stages that can be attributed to distinct growth 
laws [107j : Starting with the birth of the hole, the radius of the hole is supposed 
to grow exponentially with time, i.e. R{t) cx exp(t/r), as long as the radius R is 
smaller than R^ ~ \/hh. Experimentally, Masson et al [108] found relaxation times 
orders of magnitude larger than the largest translational reptation times expected. For 
Rc < R < R'^ ^ b, the rim is formed and viscous dissipation dominates the hole growth 
dynamics, i.e. R oc t. The special case of R ^ Rc = Vbh moreover gives an analytic 
expression for the hole growth: 

R(t) ^ ^-^(^y/h (3.22) 
T] h 

During this hnear regime, a dynamic (receding) contact angle 9d = Oy/ is formed at 
the three-phase contact line. In the subsequent stage (for hole radii approximately larger 
than the slip length b), the rim is fully established and grows in a self-similar manner. 
As discussed in section 13.2.11 this results in a characteristic R oc t^^^ growth law if large 
slip is present. That stage of self-similar growing rims is often called "mature" regime. 
These distinct regimes have been experimentally observed by Masson et al [108] and 
Damman et al [72J. Moreover, a dissipation dominated by slippage has been shown to 
be restricted to sufficiently small hole radii: If the rim has accumulated a very large 
amount of liquid (so that the height of the rim Hq is much larger than the slip length 
b, i.e. -f^o ^ b), viscous dissipation dominates again. This consequently results in a 
transition to a linear growth law for the radius with respect to the dewetting time, i.e. 
Roct. 

In this article, we mainly focus on dewetting phenomena in terms of growing holes 
or retracting straight fronts. Additionally, also the subsequent regime influenced by the 
fingering instability has attracted the interest of researchers |109j . Thereby, facets such 
as the onset of the instability and the morphology of liquid profiles have also shown 
to be sensitive to rheological properties of thin films and the hydrodynamic boundary 
condition at the solid/liquid interface |110[ lllll 11121 1113j . 

3.3.2. "Mature" holes Recently, Fetzer et al experimentally studied the dewetting 
dynamics of "mature" holes in thin PS films on two different types of hydrophobized 
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surfaces. The aforementioned assumption of linear superposition of viscous dissipation 
and slippage (see I3.2.3|l was used to gain information about the slip contribution |86j . 
Hydrophobization was achieved by the preparation of a dense, self-assembled monolayer 
of silane molecules (silanization). In that case two different silanes have been utilized, 
octadecyltrichlorosilane (OTS) and dodecyltrichlorosilane (DTS). The slip length turns 
out to depend strongly on temperature (and thereby on the melt viscosity); in fact, the 
slip length decreased with increasing temperature. Furthermore, slippage was about one 
order of magnitude larger on DTS compared to OTS. The results are in good agreement 
with the results for the slip length gained by the analysis of the shape of rim profiles of 
the corresponding holes (see 12.5.51) . The molecular mechanisms of slippage are widely 
discussed and several models have been proposed with regard to different experimental 
conditions (shear rates, polymers, surfaces, ...). This will be discussed in the next 
section. 

3.4- Molecular mechanisms of slippage 

For small flow velocities, slip lengths lower than expected have been found (see e.g. 
[114] ). This has been explained by the adsorption of polymer chains at the solid 
surface. These adsorbed chains are able to inhibit slippage due to entanglements with 
chains in the melt. At larger shear strains, disentanglement of adsorbed and melt 
chains occurs, the friction coefficient decreases rapidly and slippage is "switched on" 
[1151 11161 11171 1118] . The authors interpret the results as follows: Strongly adsorbed 
chains stay at the solid/liquid interface followed by a coil-stretch transition of these 
tethered chains. Thereby, the interaction of anchored chains (which are stretched under 
high flow velocities) with the chains in the melt are reduced and slippage is enhanced. 
This has been experimentally observed by Migler et al [114j and Hervet et al |119j . 
Thereby, the critical shear strain depends on the molecular weight of chains attached to 
the surface and on their density |39]. This situation is also called "stick-slip transition" 
and has moreover been reported for pressure-controlled capillary flow of polyethylene 
resins |120j . The authors find a characteristic molecular weight dependence of the slip 
length h oc M^'^, whereas the critical stress (Xc for the transition described before scales 
as (Tc oc Mj°-^ |121] . Also de-bonding of adsorbed chains might occur. 

Molecular dynamics (MD) studies aiming to investigate the slip length in thin, 
short-chained polymer films subject to planar shear revealed a dynamic behavior of 
the slip length 6(7) upon the shear rate according to 6 = 6*(1 — 7/7c)~^^^, where 7c 
represents a critical shear rate [122] . Slippage strongly increases as 7/7C — > 1. This 
relation was also found for simple liquids. Additionally, the authors studied the border 
case 7 — > and especially the molecular parameters affecting the slip length. 

Moreover, a special situation has been discussed in the literature: A thin polymer 
film prepared on a surface decorated with end-grafted polymer chains of the same 
species. Due to entropic reasons, an interfacial tension between identical molecules 
occurs and dewetting can take place. This phenomenon is called autophobicity or 
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autophobic dewetting |123[ 1124] . Reiter and Khanna observed that PDMS molecules 
slipped over grafted PDMS brushes. They determined a slip length on the order of 
10 /im for film thicknesses between 20 and 850 nm of the dewetting PDMS film (high 
Myj of 308kg/mol in most cases). For lower grafting densities, slippage is reduced, 
indicating a deeper interpenetration of free melt chains |123] . 

3. 5. Impact of viscoelasticity and stress relaxation 

As already introduced in 12.1.31 and furthermore theoretically described in section 12. 5[ 
viscoelasticity can be included in the description of thin film dynamics. In case of thin 
film flow with strong slippage, we already mentioned that the impact of viscoelasticity 
might be not distinguishable from slippage effects. Numerous experimental studies have 
been published by Reiter, Damman and coworkers describing the influence of viscoelastic 
effects on the flow dynamics of thin polymer films. These studies have been supported 
with theories by Raphael and coworkers. In the following, a snap-shot of recent results 
in chronological order is presented. It reflects the progress from simple to more and 
more sophisticated models and experiments. Some of the results, however, may have 
tentative character, as brought up by a very recent manuscript by Coppee et al (cf. 
section 13. 7p . 

3.5.1. Thin film rupture Reiter and de Gennes [125] pointed out that the usual 
spin-casting preparation (based on fast solvent evaporation) of thin, long-chained 
polymer films may induce a cascade of distinct states. Annealing the samples to high 
temperatures for times larger than the reptation time of the polymer (depending on 

and the temperature T) might induce a complete healing of the preparation effects. 
Non-equilibrated conformational states and residual stresses have shown to be capable 
to cause rupture of thin films |126j . Thereby, the areal density of holes appearing at a 
certain temperature above Tg has been measured: After storage at elevated temperatures 
below Tg, an exponential reduction of the hole density can be observed. Vilmin and 
Raphael were able to show that lateral stress reduces the critical value for surface 
fiuctuations initiated by an anisotropic diffusion of polymer molecules to induce the 
formation of holes |127j . 

3.5.2. Stages of distinct dewetting dynamics - Experimental results Reiter and co- 
workers furthermore studied the temporal and spatial evolution of dewetting PS fronts 
on Silicon wafers covered with a PDMS monolayer at times shorter than the relaxation 
(reptation) time r^ept in equihbrated bulk samples [126] . At the beginning of their 
experiments, the dewetted distance and the rim width w increased in a logarithmic way 
with respect to the dewetting time t, consequently V oc t~^ is found. They correlated 
a maximum of the rim width w occurring at a distinct time ti (nearly independent on 
My, and significantly shorter than x^ept) of the experiment to a change in dewetting 
dynamics, namely V oc t~^/^ (i.e. R oc \/t). 
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Even before, Damman, Baudelet and Reiter compared the growth of a hole and the 
dynamics of a dewetting straight front |72]. Thereby, strong influence of the dewetting 
geometry on the dynamics of early stages became obvious. Moreover, the authors 
correlated their findings to the shape of the corresponding rim and clearly identified 
distinct dewetting stages: At the beginning of hole growth, capillary forces dominate 
the dynamics and exponential growth is observed. In contrast to that, dewetting of 
a straight front at the beginning starts almost instantaneously at a high velocity and 
decreases as oc t^^. As already observed by Redon et al |HI], V = const and R oc t 
due to viscous dissipation is found for growing holes. For both dewetting geometries, 
Damman et al find a maximum of the rim width w which they interpret in terms of a 
transition of the shape of the rim from very asymmetric profiles towards more symmetric 
rims, while the volume still increases during dewetting. In the following dissipation- 
dominated regime, V cx t~^^^ and R oc t"^^^ (slip boundary condition) or V = const and 
R (xt (no-slip boundary condition) is found, depending on the boundary condition for 
the solid/liquid interface. Afterwards, again a constant dewetting velocity (until the 
"mature rim" regime is reached and slippage dominates) is obtained. 

Recently, Damman and Reiter determined the strain 7 from rim shapes of dewetting 
fronts, which can be written according to (12.31) by 7 = {S/h + (Tq)/G, where (Tq denotes 
residual stresses |128j . These values were one order of magnitude larger than for 
equilibrated PS films (7 = S/ {hGbuik)) and increased with increasing Mw Consequently, 
larger residual stresses (Tq or smaller elastic moduli G can be responsible for this effect. 
Indeed, the elastic modulus of a thin film is supposed to be smaller than its bulk value 
according to the fact that larger numbers for the entanglement length Mg (for larger M^) 
are expected (see section l2.1.2p . Additionally, instead of one transition time for both, 
two different transitions times t^ for the maximum of the rim width and t„ for the power 
law of the dewetting velocity (to V oc t~^/^) are found for larger molecular weights. Time 
tv is comparable to the reptation time t^ept (also for larger M^). The authors interpret 
this as follows: After ~ t^ept) the equilibrium entanglement conformation has been 
reached via interdiffusion and re-entangling of chains. Relaxation of stress, which is 
correlated to the rim shape transition and t^, occurs at shorter times. 

3. 5. 3. Stages of distinct dewetting dynamics - Theoretical models These observations 
have motivated a theoretical model by Vilmin and Raphael |129j . that takes viscoelastic 
properties of the liquid and slippage into account. The logarithmic increase of the 
dewetted distance with time can be explained if residual stress is regarded as an extra 
driving force for dewetting that has to be added to the capillary forces. Subsequent 
relaxation of the stress reduces this additional contribution. For times larger than the 
reptation time Trept, Vilmin and Raphael predict a constant dewetting velocity and 
V oc t~^^^ as soon as the "mature regime" is reached (and slippage becomes important). 
They consider a simple equation for a viscoelastic film which includes one relaxation 
time Ti (and an elastic modulus G), but two distinct viscosities: r/o is a short-time 
viscosity which describes the friction between monomers, and r/i is the melt viscosity 
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{rji <^ rjo) governed by disentanglement s of polymer chains. Theoretically, three regimes 
of distinct time- response are expected: For t < tq = tjq/G, Newtonian flow accompanied 
with a small viscosity rjo. Afterwards {tq < t < ti = rji/G), the elastic modulus governs 
the dynamics. For longer times, i.e. t > ti, Newtonian flow is again obtained, but 
with a much larger viscosity rji. During the first and the last regime, asymmetric rims 
and constant velocities are predicted (the latter velocity of course is much smaller than 
the initial velocity). The intermediate regime, governed by viscoelastic behavior, is 
expected to show V oc t^^/^ according to the fact that the width of the rim w will 
increase proportional to the dewetted distance D, i.e. w oc D (in contrast to the square- 
root dependence in case of viscous flow). Including residual stresses into their model, 
Vilmin and Raphael obtained the aforementioned experimentally observed V oc t"^ law 
instead of oc t~^^'^. Recently, Yang and co-workers experimentally quantified the 
molecular recoiling force stemming from non-equilibrium chain conformations. They 
obtained values that are at least comparable to (or even larger than) the dispersive 
driving forces [130] . 

3.5.4- Further remarks Concerning studies and models including viscoelastic effects, 
one has to bear in mind that stress relaxation dynamics and relaxation times in thin films 
have been shown to be different from bulk polymer reptation values. As mentioned in 
section 12.1.21 and indicated by the characteristic relaxation times ti of elastic constraints 
being significantly shorter than bulk values, the entanglement length has shown to be 
significantly increased in case of thin films. 

Recently, Reiter and coworkers extended their studies with regard to viscoelastic 
dewetting on soft, deformable substrates [100] . The essential result was that transient 
residual stresses can cause large elastic deformations in the substrate which almost stop 
dewetting for times shorter than the relaxation time Trept of the polymer film. For times 
longer than r^ept, the elastic behavior and the elastic trench in the deformable substrate 
vanishes. 

Vilmin and Raphael applied their model for viscoelastic liquids and residual stresses 
to the hole growth geometry [131] with regard to the early stage (exponential growth 
for Newtonian liquids). They discovered a very fast opening regime followed by a slow 
exponential growth of the radius of the holes. 

3.6. Non-linear friction 

Up to now, concerning energy dissipation at the solid/liquid interface, exclusively linear 
friction has been considered. As described in section 13.41 often not only smooth and 
passive (non-adsorbing) surfaces have been experimentally considered, but also grafted 
or adsorbed polymer layers are used as a support for dewetting experiments. These 
systems motivate the theoretical treatment of non-linear friction, as described in the 
following subsection. 
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3.6.1. Theoretical model To cover non-linear friction, Vilmin and Raphael introduced 
a friction force per unit area that is linear below a certain transition velocity Vr and 
non- linear above |131j : 

Fs = iVr{v/Vrf-'', V > Vr, (3.23) 

where r denotes a so-called friction exponent. 

Consequently, the effective (velocity-dependent) friction coefficient ^e// can be written 
in terms of 

^effiv) = aVr/vY (3.24) 

and a velocity-dependent slip length b{v) can be defined as 

Kv)=r^/^eff{v). (3.25) 

An increasing velocity leads to a decreasing friction coefficient and to pronounced 
slippage. Thereby, the friction at the solid/liquid interface influences the power law 
of the velocity decrease. Assuming non-linear friction in the intermediate regime (which 
is governed by the viscoelastic behavior) gives (according to w oc D) V oc t~^^^'^~^\ 
Including residual stress o"o to this model leads to a formula for the maximum width 
of the rim (depending on r and (Jq). Experimentally, values for r between (for low 
My^) and 1 (for large M^) have been found |128j . The dependence of the nonlinearity of 
friction upon molecular weight My^ could be explained by the influence of chain length 
on slippage (see f l3.12p ). 

3.6.2. Variation of substrate properties Hamieh et al focused on the frictional behavior 
of dewetting viscoelastic PS films on PDMS-coated (irreversibly adsorbed) silicon wafers 
|132] . Thereby, thickness (via the PDMS chain-length) and preparation of the PDMS 
support (via the annealing temperature) were varied. In summary, the observations 
are consistent with the aforementioned (section 13. 5p experiments by Reiter, showing a 
characteristic time ti for the change in dewetting dynamics and rim shape (from highly 
asymmetric towards a more symmetric equilibrated shape). Again, this transition is 
interpreted by the Laplace pressure that overcomes elastic effects. Probing the maximum 
rim width enables the authors to identify the impact of the friction coefficient or the 
friction exponent r: If they prepare thicker PDMS layers, the result is a larger maximum 
rim width. This result can be explained in terms of a small increase of r. Consequently, 
the velocity-dependent slip length h{v) increases (see (13.251) ): Thicker PDMS-layers 
lead to more slippage. Concerning the characteristic time for stress relaxation ti, no 
significant influence of the preparation procedure (annealing temperature of the PDMS 
coating) and thickness of the PDMS layer have been found. 
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3. 7. Temporal evolution of the slip length 

Most of the aforementioned dewetting experiments of thin (viscoelastic) PS films are 
based on supports consisting of a PDMS layer prepared onto a Si wafer. These PDMS 
surfaces, which were assumed to be impenetrable by PS chains, turned out to be 
less ideal. Recently, the aforementioned fast decay of the dewetting velocity and the 
maximum of the rim width upon dewetting time (ascribed to relaxation of residual 
stresses, see section 13751) . has also been observed in case of low molecular weight PS [133] . 
Furthermore, neutron reflectometry experiments revealed an interdiffusion of polymer 
chains at the PS-PDMS interface below the PS bulk glass transition temperature. The 
dewetting velocity accelerated as the brush thickness is increased. In view of energy 
dissipation due to brush deformation (pronounced by increasing thickness), this is 
contradictory to the slower dewetting velocities expected in that case. 

Ziebert and Raphael recently investigated the temporal evolution of the energy 
balance (viscous dissipation and sliding friction) of thin film dewetting by numerical 
treatment, especially concerning non-linear friction [L3_4] . They point out that both 
mechanisms have different time dependencies and propose that simple scaling arguments 
such as the mass conservation w (x D should be revisited. In case of non-linear friction, 
viscous dissipation is even more important than sliding friction for times larger than ri, 
whereas for linear friction both mechanisms are approx. equally important. The stage 
of "mature" rims afterwards is again dominated by friction at the solid/liquid interface. 
Using scaling arguments and numerical solutions of the thin film model for viscoelastic 
liquids |129j . they moreover showed that the aforementioned dewetting characteristics 
(fast decay of the dewetting velocity, maximum of the rim width in course of dewetting 
time) can be explained by the temporal decrease of the slip length during the experiment 
instead of stress relaxation |135j . Therefore, roughening of the PS-PDMS interface (as 
detected by neutron reflectometry) or potentially also an attachment of very few PS 
chains to the silicon substrate might be responsible, if very low PDMS grafting densities 
are prepared. 

The attachment of melt chains to the substrate has been recently considered by 
Reiter et al |136] . They pointed out that the driving force is reduced by a certain pull- 
out force for molecules attached to the surface as they get stretched while resisting to be 
pulled out. Consequently, this force (per unit length of the contact line) can be written 
in terms of Fp = vLf* , where v represents the number of surface-connected molecules, 
L the length according to the stretching and /* the pull-out force per chain. 

4. Conclusions and outlook 

To conclude, dewetting experiments can be regarded as a very powerful tool to probe 
rheological thin film properties and frictional mechanisms. The validity of the no-slip 
boundary condition at the solid/liquid interface, for a long period of time accepted as a 
standard approach in fluid dynamics, fails. Even for Newtonian liquids such as polymers 
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below their critical length for entanglements substantial amount of slippage has been 
found. For larger molecular weights, slippage can be even more pronounced if chain 
entanglements come into play. 

We have reviewed in this article theories and experiments characterizing the statics 
and dynamics in thin liquid polymer films, focussing on energy dissipation mechanisms 
occurring during dewetting. Experiments are relatively easy to perform, since usually 
no very low or very high temperatures are needed, neither are high-speed cameras 
necessary. However, careful preparation, preferably in a clean-room environment, of 
thin films is inevitable. Concerning the experimental data, diligent interpretation and 
consideration of all relevant parameters are of essential importance: Molecular weight, 
film thickness (also with regard to the molecular dimensions of a polymer coil in the 
respective melt), dewetting temperature, melt viscosity, dewetting velocity and capillary 
number, residual stresses, relaxation times and ageing time represent parameters that 
are sometimes coupled and not easy to disentangle. Their unique impact on the friction 
coefficient or the slip length on a specific, ideally non-adsorbing and non-penetrable 
substrate (to reduce the number of parameters of the system) is not always easy to 
identify. In particular, the specific stage of dewetting (early or mature stage) and the 
dewetting geometry (holes or straight fronts) represent further facets and playgrounds 
for experimentalists and theoreticians that, under careful consideration, enable to gain 
insight into rheological or frictional mechanisms. The dynamics and the morphology 
of the fingering instability can provide additional access to the solid/liquid boundary 
condition and the rheology. 

To obtain a deeper understanding of the molecular mechanisms at the solid/liquid 
interface is one of the main tasks in micro- and nanofiuidics. In this context, we like to 
highlight two essential questions and possible pathways to answer them: 

a) The liquid: What is the impact of the polydispersity of the liquid on dewetting? 
Dewetting studies of polymer melts by adding a second chemical component have shown 
to infiuence slippage |137j . However, also the driving force is changed due to the 
difference in chemical composition for different species. To overcome this problem, 
the influence of chain length distribution on dewetting can be probed by studying 
polymer mixtures instead of monodisperse polymer melts [138] . Concerning a theoretical 
approach to this question, dissipative particle dynamics (DPD) simulations enable to 
locate energy loss in the rim and reveal interesting results in case of two immiscible fluids 
of different viscosity: In case of a low- viscosity layer at the solid/liquid interface, faster 
dewetting dynamics is found that is attributed to a lubrication effect, i.e. the sliding of 
the upper high- viscosity layer |139j . Finally, these aspects lead to the fundamental 
discussion, whether "apparent" slip, possibly induced by the formation of a short- 
chained layer of low viscosity, is present. 

b) The substrate: What is the impact of the molecular structure of the substrate? 
As indicated before, the set of parameters concerning the topographical and/or chemical 
structure of the support is large. Besides parameters such as surface roughness and 
surface energy, experiments can be performed on substrates e.g. decorated with 
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an amorphous coating, a self-assembled monolayer or ever grafted polymer brushes 
of the same or different species. Scattering techniques (using neutrons or X-rays) 
provide access to the solid/liquid interface, complementing dewetting experiments and 
confirming proposed mechanisms of slippage. Simulations based on molecular dynamics 
(MD) of near-surface flows can help to compare experimental results from dewetting 
studies to molecular parameters, easily tunable in theoretical models, and structural 
changes [140[ I141j . Moreover, the evolution of coarse-grained polymer brush/melt 
interfaces under flow has also been identified as a potential application of MD studies 
[142[ I143j . In the end, all these facets help to obtain a universal picture of sliding 
friction which can potentially lead to surfaces precisely tailored for special microfluidic 
applications. 
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